Synchronization of Random Linear Maps 
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We study synchronization of random one-dimensional linear maps for which the Lyapunov exponent can be 
calculated exactly. Certain aspects of the dynamics of these maps are explained using their relation with a ran- 
dom walk. We confirm that the Lyapunov exponent changes sign at the complete synchronization transition. 
We also consider partial synchronization of nonidentical systems. It turns out that the way partial synchroniza- 
tion manifests depends on the type of differences (in Lyapunov exponent or in contraction points) between the 
systems. The crossover from partial synchronization to complete synchronization is also examined. 

PACS numbers: 05.45.Xt 



I. INTRODUCTION 

Synchronization of chaotic systems is a subject of current 
intensive study d. To a large extent this is due to its var- 
ious applications, ranging from laser dynamics |2] to elec- 
tronic circuits 01, chemical and biological systems 0], se- 
cure communications |5], etc. But there is also a pure theo- 
retical interest in this phenomenon, that is related perhaps to 
its counterintuitive nature: how it is possible that chaotic (i.e., 
by definition unpredictable) systems can be synchronized and 
thus brought under some 'control'. And, even more puzzling, 
noise can play the role of the synchronizing factor. Indeed, 
early reports |6] that sufficiently strong noise can completely 
synchronize two identical chaotic systems were initially met 
with scepticism and attributed to finite precision of computa- 
tions 1 7] or to biased noise [8]. However, more recent exam- 
ples show this effect even for unbiased noise |9]. 

Since real systems are typically nonidentical, complete syn- 
chronization is difficult to achieve. It is an interesting problem 
to examine whether noise can induce some sort of 'weaker 
synchronization' in nonidentical but relatively similar sys- 
tems. Recent works do show the existence of such partial 
synchronization floHTlll . 

One important problem of the theoretical and numerical 
studies of synchronization is how to detect it. This problem 
is essentially solved for the complete synchronization of two 
identical systems that are described by variables x and x', re- 
spectively. In this case, for the synchronized state the dif- 
ference \x — x'\ equals zero, while it remains positive in the 
unsynchronized state. Moreover, the transition between these 
two states is accompanied by the change of the sign of the 
largest Lyapunov exponent, that becomes negative in the syn- 
chronized state. 

However, for the partial synchronization this problem is 
much more subtle. In this case, the two systems are not identi- 
cal and the difference \x — x'\ always remains positive. It has 
already been noted for some models with continuous dyna- 
mics that partial synchronization manifests through changes 
in the probability distribution of the 'phase difference' II ill . 
In addition to that, in the partially synchronized phase the so- 
called zero Lyapunov exponent becomes negative [ 1(1 Hill . 

Studies of synchronization rely, to a large extent, on nu- 



merical calculations. Precise estimations of Lyapunov expo- 
nents or probability distributions (invariant measures) consti- 
tute very often demanding computational problems. To fur- 
ther test the already accumulated knowledge on synchroniza- 
tion, it would be desirable to find models where at least some 
of these properties could be computed analytically. 

In the present paper we examine synchronization of ran- 
dom one-dimensional linear maps. For such maps one can 
easily find the exact Lyapunov exponent and locate the point 
where it changes sign. Numerical calculations for two iden- 
tical systems confirm that this is also the point where a com- 
plete synchronization transition takes place. We briefly report 
on a correspondence between such maps and a random walk 
process, that allows for a simple interpretation of the initial 
stages of the evolution of the maps. 

We also examine the partial synchronization of noniden- 
tical maps. It is seen that the way partial synchronization 
manifests depends on the type of difference between the non- 
identical systems. In a certain case, the difference in the lo- 
cation of contraction points of the maps is imprinted in the 
probability distribution at the partial synchronization transi- 
tion. When the difference between two subsystems S tends to 
vanish, partial synchronization approaches the complete syn- 
chronization. Due to the exact knowledge of the complete 
synchronization point, one can examine some details of this 
crossover. In particular, it is shown that for 5^0, vanishing 
of the synchronization error is very slow ~ (— l/log 10 S). 



II. RANDOM 2-MAPS 

First, let us consider the simplest example of a random lin- 
ear map 



c„+i = fi(x n ), i = 0, 1 and n = 0, 1, 



(1) 



where at each time step n one of the maps fo or /i is applied 
with a probability p and (1 — p), respectively. The maps are 
defined as fo(x) = ax mod(l), and fi(x) — bx, where 
< x < 1 and a > 1, < b < 1. Related models 
have already been examined in the context of on-off intermit- 
tency rjUl. advection of particles by chaotic flows 0, 
and others I151 fl6ll . Some aspects of synchronization were 



also studied for piecewise linear random maps, but both the 
Lyapunov exponent and the location of the synchronization 
transition were determined only numerically 11711 . 

For the map Q it is elementary to calculate exactly its Lya- 
punov exponent A defined as 



1 

li m TT / logi 

71 = 



dx 



n+l 



dXr 



(2) 



Indeed, since both /o and f\ have constant derivatives, one 
immediately obtains 



A = P logio a+(l-p) log 10 b 



It follows that A changes sign at 



P = Pc 



-logio b 
log 10 (a/b) 



(3) 



(4) 



It is easy to understand this result. For p > p c the expand- 
ing (chaotic) map fo prevails over the contracting map /i and 
the overall behaviour is chaotic with A > 0. The opposite 
situation takes place for p < p c and the map Q contracts to 
x = 0. At p = p c and in its close vicinity the map Q exhibits 
intermittent bursts of activity Jl2lfT3ll . but we will not focus 
on such a behaviour in the present paper. 

To study synchronization one can make two runs {x„ } and 
{x' n } of iterations of map ([0 starting each time from (slightly) 
different initial conditions xq and x' but with the same real- 
ization of noise, i.e., with the same sequence of maps /j. Then 
one measures the synchronization error w n defined as 



■'' u I) • 



(5) 



where (...) represents the mean over the realizations of the 
noise. Moreover, we introduce the steady-state average w = 
lim n _ >00 , i<; n . In the synchronized state w = while it is pos- 
itive in the unsynchronized state. Our numerical results for 
a = 1/6 = 3/2 (not presented here) show that w vanishes at 
p = 1/2, i.e., the point where the Lyapunov exponent changes 
sign. 

However, due to the fact that it contracts to x = for p < 
p c , the map Q is not quite suitable to study synchronization 
of chaotic systems. For such a purpose it would be desirable 
to have a map with a somehow more complex behaviour in the 
regime with a negative Lyapunov exponent. 

Nevertheless, the simplicity of the map Q allows us to get 
some additional insight into its dynamics. Let us fix the initial 
points of our maps as, e.g., xq < x' <C 1. As we shall see, 
such a choice results in a certain transient regime that can be 
understood using an analogy with a random walk process. 

Numerical evaluation of w n shows that it exhibits three 
types of behavior as a function of time n (Fig. 0. 

(a) For p > p c (positive Lyapunov exponent) after an ini- 
tial exponential increase, w n saturates and acquires a constant 
nonzero value. 

(b) For p < p c , w n decreases after an initial exponential in- 
crease. The asymptotic exponential decay is consistent with 
the (negative) Lyapunov exponent (|3}. 
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FIG. 1: The logarithm of the synchronization error w n as a func- 
tion of time n for the map (1} with a = 1/6 = 3/2, for different 
values of p. The mean in w n was taken over 10 10 realizations of the 
random trajectories, and the initial points were xo = x' /2 — 10~ 9 . 
The diamonds locate the maximum rp of the corresponding proba- 
bility P(n; p, yo — 50). The inset shows a typical such first-passage 
probability P(n;p = 0.55, yo = 50) for the corresponding random 
walker to go from yo to in n steps. 



(c) For even smaller values of p the synchronization error w n 
decreases exponentially already from the beginning. 

To understand the initial behaviour of w n , let us recall that 
the iteration of x n and x' n starts from very small values and 
for a certain number of iterations the mod(l) part of the map 
Q does not play any role. Consequently, one has 



where 



w n = (a) n w , 



(a) = pa + (1 — p) b. 



(6) 



(7) 



Correspondingly, the initial decrease or increase of w n de- 
pends whether p is smaller or greater than (1 — b)/ (a — b). 

After a certain time, the mod(l) part of the map comes 
into play and the initial behaviour of w n is replaced by a 
different one. Namely, for p > p c , w n saturates at a nonzero 
value and it decays exponentially for p < p c . We shall return 
to this point at the end of this section. 

To estimate the time scale r when the inital behaviour (|SJl 
changes, we relate our map to a random walk process. For 
simplicity, we consider the case b — 1/a. In this case there 
is a one-to-one correspondence between the stochastic vari- 
able x n and y n — ~log a (x n ). Multiplication of x n by a or 
1/a corresponds to the decrease or increase of y n by unity. 
Consequently, y n is nothing else but the position of a random 
walker on a lattice of unit spacing with transition probabili- 
ties p to the left and (1 — p) to the right. The correspondence 
with the random walk holds as long as the mod( 1 ) part of the 
map is not applied, i.e., the walker does not cross log a (l) = 0. 
The above random-walker problem has two characteristic time 
scales connected to the first-passage process 1 18] from its ini- 
tial position to 0, and one can expect that one of them is re- 
lated to t [which is, recall, the time scale on which the initial 
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behavior of the map (jfji is altered by the mod(l) part of the 
map], (i) First, there is the mean first passage time tm- But 
for p < p c , tm is infinite - contray to r, and therefore one 
cannot use tm as a measure for r. (ii) Secondly, there is the 
time moment Tp when the probability that the walker hits 
for the first time is maximal. The probability distribution of 
the first hit is known to be ifTsIl : 



P(n;p,yo) = — 

n 



n-y p 2 (1-p) 2 



(8) 



where yo is the initial position of the walker and n is the 
number of steps; the binomial coefficient is to be interpreted 
as zero if (n — yo)/2 is not an integer in the interval [0, n]. 
Our numerical simulations show (see Fig.Q that the value of 
n = Tp for which P(n]p,yo) in Eq. attains a maximum 
offers a reasonable estimation of r. One cannot expect to get 
a more precise estimate of r, since the change of the initial 
behavior (|6} of the map towards the asymptotic one is a grad- 
ual process, that involves all the trajectories of the equivalent 
random walker that hit (for the first time) before as well as 
after Tp. 

Let us notice that the bias of the random walk is related 
with the sign of the Lyapunov exponent and therefore with the 
asymptotic behaviour of our map. Indeed, forp < p c = 1/2 
the random walk is biased toward +oo, that translates into an 
exponential decay of w n . Our numerical simulations for the 
longer time regime suggest that this decay is governed by the 
Lyapunov exponent A 0. Recall that A is known to govern 
the evolution of the typical value of \x n — x' n \ (see, e.g., lllrjlo . 
Thus in the long-time regime the mean w n and the typical 
value of \x n — x' n \ behave identically. On the other hand, 
they are clearly different in early-time regime. It means that 
in this regime w n is strongly influenced by rare events, i.e., 
unlikely excursions of the random walker against the bias. For 
p > p c ~ 1/2 the random walk is biased toward and, since 
the map is bounded, w n saturates at a positive value. 



HI. RANDOM 3-MAPS 

As already mentioned, the map Q has a trivial behaviour 
for p < p c and is not suitable to study synchronization of 
chaotic systems. In this context the following 3-map version 
is more interesting: 



x n +i = fi(x n ), i = 0,1,2, 



where /o(a;) = ax mod(l), fi{x) = bx, and f2(x) = bx 
(1 — b). The maps /o, fx, and f 2 are applied at random with 
probabilities p, (1 — p)/2, and (1 — j>)/2, respectively. It is 
easy to show that for such a random map Eqs. still 
hold. 

Numerical evaluation of the synchronization error w — 
lim ?woo w„ for the map (|9} with a = 1/b = 3/2, based on 
N = 10 8 iterations, is shown in Fig.|2] Forp > p c = 1/2 the 
Lyapunov exponent A is positive and the system is not syn- 
chronized (w > 0). For p < 1/2 we have A < and the 
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FIG. 2: The synchronization error w as a function of p calculated for 
the map |9j with a = 1/6 = 3/2. The unnormalized probability 
distribution P(x) of visiting a point x for p = 0.4 (synchronized 
phase) is shown in the inset. The data were collected in 10 4 bins. 
Complex probability distributions are known to appear also for some 
other random linear maps 1 1 311 - 



system synchronizes (w — 0). But this time the behaviour for 
p < 1/2 is much more complex. The two maps /i and fa 
are still contracting ones, but to two different points (0 and 1). 
Since they are applied randomly, the system, while remaining 
synchronized, irregularly wanders throughout the whole in- 
terval (0, 1). The probability density Fix) of visiting a given 
point x, that is shown in the inset of Fig.|2j confirms such a 
behaviour. 

Random linear maps can be also used to study partial syn- 
chronization that might occur for two dynamical systems that 
are not identical, although relatively similar. Generally, one 
considers the following pair of maps 



x n+l — fi{, x n)-i x n +l — fi i.Xn)i ^ — 0,1,2, 



(10) 



for which {/o, /o}, {fi, fi}, and {/ 2 , J' 2 } are applied 
at random with probabilities p, (1 — p)/2, and (1 — p)/2, 
respectively. To complete the definition we have to specify 
the functions fi and i = 0, 1,2. We present below the 
results for two particular choices. These ones correspond, 
respectively, to a perturbation of the Lyapunov spectrum - 
case (A) - and of the attractor - case (B). 



(9) Case (A): We choose 



+ ffa) = a' x mod(l), f[(x)=b'x, &{x) = b' x + (1 - b') 

(11) 

where a' = 3/2 + 5, b' = I /a' and S = 1CT 3 . Functions / , 
fx, and /2 are defined as for the map Q with a = 1/b = 3/2. 
For such a choice, the Lyapunov exponent of the map in 
Eq. il It is a linear function of p that also changes sign at 
p = p c = 1/2 albeit with a different slope. Assuming that a 
partial synchronization transition is also governed by the Lya- 
punov exponent, we might expect that such a transition, if it 
exists, takes place at p = 1/2. For the pair of maps i ll Oi l. 
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FIG. 3: The synchronization error w and the probability s that \x — 
x'\ < 10 -4 as a function of p for map 

{10}-{ID- The unnormalized 
probability distribution P(\x — x'\) is shown in the inset. 
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FIG. 4: The synchronization error w and the probability s that \x — 
x'\ < 10 -4 as a function of p for map i 101 - i 1 2i . The unnormalized 
probability distribution P(\x — x'\) is shown in the inset. 



w n is defined as in Eq. (0 but x n and x' n evolve according 
to Eq. JlOb . Numerical calculation of w = lim„_ +00 u' rl for 
N = 10 s iterations shows (Fig. [3} that it remains positive 
and smooth around this value. However, there is a more sub- 
tle change in the system at or around p = 1/2. Indeed, the 
probability distribution of P(\x — x'\) shows a pronounced 
peak in the p < 1/2 domain (inset of Fig. This is yet 
another example that shows that partial synchronization man- 
ifests through a change in the probability distribution |11]. 
However, P(\x—x'\) was calculated forp = 0.4, 0.5, and 0.6, 
i.e., values that are relatively far from each other. For smaller 
differences between the values of p the corresponding curves 
are getting similar and it is not clear to us whether there is a 
qualitative change in the probability distribution P(\x — x'\) 
that could locate precisely the partial synchronization transi- 
tion. If not, it could mean that in this case either there is not a 
well-defined partial synchronization transition, or we are not 
looking at the right quantity to detect it. 

For further comparison, we also calculated the probability s 
that the difference \x — x' | remains smaller (if initially so) than 
a given value e = 10~ 4 during the iterations of the system 
J10I . One finds that s is a monotonous function of p (Fig.[3j. 

Case (B): We choose 

f' Q {x) = a'x mod(l), f[{x) = b' x, 

ti(x) = b'x+(l-b')-b" , (12) 

where a' = 3/2, b' = 2/3 and b" = 10~ 3 . Functions / , 
/i, and /2 are defined as in case (A). For such a choice the 
Lyapunov exponents of both systems are the same and they 
change sign at p = p c = 1/2. The only difference is that the 
function /^(s) has a contracting point at 1, while /^(a/) at a 
slightly smaller value x' — 0.997. 

In this case the synchronization error w behaves similarly 
to case (A) (Fig.|4}- However, the probability distribution of 
P ( | x — x' | ) has a little bit different shape. In particular, even at 



1.14 




1.08 1 1 1 1 

0.0025 0.003 0.0035 0.004 

lx-x'l 



FIG. 5: The unnormalized probability distribution P(\x — x'\) cal- 
culated for map 1 101 - 1 121 at p — 0.5. The maximum at \x — x'\ = 
0.003 is due to the shift in the location of contracting points. Off the 
synchronization transition (p = 0.5) the location of the maximum 
changes. 



p = 1/2 it has a certain peak which increases with decreasing 
p. The maximum of the peak at p = 1/2 is most likely located 
&t\x-x'\ = 0.003 (Fig.|3 and this is related with the dif- 
ference of the contracting points of the functions / 2 (x = 1) 
and (x' = 0.997). We checked that off the p = 1/2 point 
the maximum shifts away from \x — x'\ = 0.003. Moreover, 
the probability s that the difference \x — x'\ is smaller than 
e = 10 -4 shows a maximum at p = 1/2 (Fig.|4j. 

Comparison of partial synchronization in cases (A) and (B) 
revealed important differences between them. In case (A), 
where we perturbed the Lyapunov exponent, the probability 
distribution P(\x~x'\) develops a peak around/at \x — x'\ = 
in the partially synchronized state that presumably exists for 
p < 1/2. Such a feature is similar to those already reported in 
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FIG. 6: The synchronization error w as a function of the mismatch 
S for maps 11 W . The inset shows the data only at the partial synchro- 
nization transition p = 0.5. 

the literature full . However, it is not clear to us whether in this 
case such a change can be characterized more quantitatively 
so that a well-defined transition point exists. A different be- 
haviour takes place in the case (B) with perturbed contracting 
points. In this case the peak of the probability distribution is 
shifted by a value that at the partial synchronization transition 
p = 1/2 (as deduced from the vanishing of the Lyapunov ex- 
ponent) curiously matches the shift of the contracting points. 
In addition to that, the probability of being in a state with a 
small difference \x — x'\ has a maximum at p = 1/2 - in 
drastic contrast with case (A). 

For simple maps as those examined in the present paper one 
has a complete knowledge of the Lyapunov exponents and 
fixed points. This is usually not the case for more compli- 
cated dynamical systems like Lorentz or Rossler equations. 
Perturbing some parameters of these equations one usually 
modifies both the Lyapunov spectrum and the attractors. It is 
henceforth possible that partial synchronization in such sys- 



tems combines some features of both cases (A) and (B) we 
discussed. 

When the difference between nonidentical systems van- 
ishes, partial synchronization is replaced by complete syn- 
chronization. We shall now briefly examine such a situation. 
In particular we consider two nonidentical systems as those 
described in the case (A) above, with varying difference 5. 
We calculated the synchronization error w as a function of 8 
for p — 0.4 and p = 1/2. In both cases we expect that w 
vanishes for 8 — > 0. Figure |6] confirms such an expectation, 
but it also reveals that for p = 1/2 the convergence to zero 
is much slower than for p — 0.4. While for p — 0.4 the syn- 
chronization error seems to vanish as w ~ 8, the inset in Fig. [6] 
suggests that at the synchronization transition the vanishing is 
most likely logarithmic w ~ (— l/log 10 8). 

IV. CONCLUSIONS 

In the present paper we examined synchronization of one- 
dimensional random maps. Our results confirm that a com- 
plete synchronization transition coincides with the change of 
sign of the Lyapunov exponent. We also showed that the way 
the partial synchronization manifests depends on the type of 
difference between the two nonidentical systems. It would 
be desirable to explain the origin of the logarithmically slow 
crossover of the partial synchronization to the complete syn- 
chronization. It would be also interesting to explain why the 
location of the maximum of the probability distribution at the 
partial synchronization transition discussed in case (B) above 
matches exactly the shift of the contracting points. 
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